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1.0  INTRODUCTION 


The  long  range  objectives  of  this  research  are  to  identify  and  evaluate  the  heat  transfer  char¬ 
acteristics  of  evaporating  thin  liquid  films. 

The  optimal  design  of  a  heat  exchanger  is  the  one  that  provides  the  highest  heat  transfer  rate 
at  the  lowest  overall  cost.  Deijaguin  [1]  introduced  the  concept  of  the  “disjoining  pressure”  to 
describe  the  effect  of  Van  der  Waals  dispersion  forces  in  very  thin  films  and  their  role  in  the  dry¬ 
out  of  liquid  films.  Below  a  critical  film  thickness,  a  non-evaporating  absorbed  layer  exists.  As 
the  liquid  film  increases  in  thickness,  the  disjoining  pressure  gradient  enhances  flow  into  the  film 
which  they  suggested  could  increase  the  evaporation  rate  several  fold.  Potash  and  Wayner  [2] 
experimentally  determined  two  distinct  regions  of  the  evaporating  extended  meniscus  on  a  verti¬ 
cal  flat  plate.  Miller  [3]  studied  the  stability  of  liquid-vapor  interfaces  moving  as  a  result  of  phase 
transformation  or  mass  transfer.  Das  Gupta  et  al  [4]  used  the  Kelvin-Clapeyron  change  of  phase 
heat  transfer  model  to  evaluate  experimental  data  for  an  evaporating  meniscus. 

In  the  present  study,  attention  is  focused  on  the  analysis  of  the  following  aspects  of  thin  liq¬ 
uid  films: 

1.  The  effect  of  an  electric  field  on  the  stability  of  a  non-evaporating  thin  film. 

2.  The  effect  an  of  electric  field  on  evaporating  thin  films. 

The  formulation  of  some  of  these  problems  will  be  summarized  in  this  report 
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2  »  ^“‘^^OSTATIC  FIELD  ON  NON-EVAPORATING 

2  liquid  film  as  it  flows  down  an  inclined  plane  was  studied  in  references 
[5-10].  AU  these  research  workers  ignored  the  London  Van  der  Waals  attractions  between  the 

solid  surface  and  the  liquid  and  therefore  their  results  cannot  be  directly  applied  to  the  case  of  thin 
hquid  films. 

Our  aim  here  is  to  address  the  question  of  how  the  thin  liquid  film  and  an  electrostatic  field 
interact  We  are  interested  in  the  specific  working  regimes  of  the  parameters,  where  it  will  be  pos¬ 
sible  to  stop  a  rupture  (or  dryout)  of  the  film.  This  will  be  done  by  solving  the  equations  of  film 

motion  in  the  presence  of  an  electric  field.  Questions  concerning  the  film  stability  will  be 
addressed. 

We  consider  the  flow  of  a  thin  liquid  film  down  an  inclined  plane  under  gravity.  The  plane 
is  assumed  to  make  an  angle  P  with  the  horizontal.  We  choose  x  and  y  directions  to  be  parallel 
^d  normal  to  the  plane,  respectively  as  shown  in  Figure  1.  We  assume  that  the  characteristic 
thickness  of  the  film  to  be  d  and  the  length  scale  parallel  to  the  film  to  be  L.  The  aspect  ratio 
d/L .  For  a  thin  film,  ^  «  1 .  The  distance  from  the  charged  foil  and  the  plane  is  H. 


Figure  i.  Flow  Model  for  the  Thin  Film  Flow 
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The  electric  field  is  determined  by  solving  Laplace’s  equation. 

V2(|)=  0  (1) 

where  (|>  (x,  y)  is  the  electric  potential.  The  boundary  conditions  are 

(|)UiD  =  <D(.r);(t>U0)  =  0  (2) 

Along  y  =  hix,  t)  we  have  the  boundary  conditions  that  the  tangential  electric  field  and  the  nor¬ 
mal  displacement  field  are  continuous.  It  may  be  noted  that  y  =  /i  {x,  t)  is  unknown,  so  that 
solution  of  the  electrostatic  problem  is  coupled  to  the  dynamics  of  the  film. 

The  liquid  film  is  governed  by  the  Navier-Stokes  equations.  The  Uquid  layer  is  assumed 
thin  enough  that  Van  der  Waals  forces  are  effective  and  thick  enough  that  a  continuum  theory  of 
the  liquid  is  ^plicable. 

We  assume  that  the  liquid  is  incompressible.  The  governing  equations  and  boundary  condi¬ 
tions  are  made  dimensionless  by  using  the  following  scales:  length  in  y-direction  ~d,vD.  ;c-direc- 
tion  »  L ,  velocity  in  j:-direction  =  ,  velocity  in  y-direction  unit  of  time  =  L/U^ ,  unit 

ofpressurepC/2, and  unitofelectric  field  =F.  We  take  p  as  the  fluid  density,  thedielectric 
constant  and  p.  the  viscosity.  The  continuity  equation  becomes 


du  9v 
dx^  dy  ~ 


The  momentum  equation  becomes 

dif  du\  dp  1  9w 

In  the  above  equations,  u  and  v  are  the  velocity  components  in  a:  and  y  directions  respectively,  p  is 
the  pressure  and  \|/  is  the  dimensionless  potential  function  representing  the  van  der  Waal’s  forces. 
We  follow  Williams  and  Davis  [10]  and  write  a  modified  expression  for  \|/: 

V  =  Ah-^  (5) 

where  A  is  related  to  the  Hamaker  constant  A'  as  A  =  .  We  have  introduced  the  Reynolds 

number.  Re  =  pU^d/iL  and  the  Froude  number,  Fr  =  UQ/Jgd. 

The  boundary  conditions  along  the  solid  plane  wall  are  given  by: 

y  =  0;M  =  v  =  0  n\ 
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At  the  fluid  interface,  we  have  the  kinematic  condition: 


y  =  hix,t): 


dh  dh 


dt 


+  =  V 


The  continuity  of  tangential  stress  on  the  interface  requires 

The  continuity  of  normal  stress  at  the  interface  y  =  h{x,t)  becomes 


(8) 


(9) 


(10) 


^  =  i6iirz?o 

=  dielectric  constant  of  the  fluid,  dielectric  constant  for  vapor  assumed  =  1 
Eq  =  electrical  permittivity  of  free  space 

,  =  normal  and  tangential  components  of  the  electric  field  in  the  vapor  at  the  inter¬ 
face 

Ca  =  is  the  capillary  number 

The  left  hand  side  of  Eq.  (10)  includes  the  term  Ah~^  to  describe  disjoining  pressure  effects 
explicitly  in  the  interfacial  boundary  condition.  There  is  still  some  question  about  whether  this 
should  be  included.  Williams  and  Davis  [10]  did  not  include  this  term  in  a  similar  analysis.  This 
does  not  change  the  leading  order  values  for  w;,  vj,pj,  etc.  Equations  (3)  -  (10)  determine  the 

motion  of  the  liquid  film.  Our  aim  here  is  to  solve  for  the  stability  of  the  liquid  film  while  includ¬ 
ing  the  effect  of  Van  der  Waais  forces  and  an  applied  electric  field. 

We  now  apply  the  long-wave  theory  to  study  the  stability  problem.  When  the  layer  is  thin¬ 
ner  than  a  critical  value,  small  disturbances  begin  to  grow.  These  waves  have  wavelengths  much 
larger  than  the  mean  thickness  of  the  layer.  Defining  a  small  parameter  k  that  is  related  to  wave 
number  of  such  disturbances,  we  may  rescale  the  governing  equations: 

X  ^  Kx;  Y  =  y;x  =  Kt  (H) 
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we  now  assume  the  following  expansions  for  the  flow  field: 

U  =  Uq+  KMj  +  +  0  (k^) 

V  =  k[V(,+  KVj+k2v2  +  0(k3)  ] 

1 

P  =  +  ^  Pi  +  0 (k3)  ]  (12) 

\|;Q  =  K\j;  =  0(l)  asK-»0 

Neglecting  the  variation  of  p  with  y  and  substituting  expressions  (12)  into  equations  (1)  -  (10),  we 
get 


^0  = 


-Re 


“0  =  [-|^sinP  +  Re 

^■'^9X2  >16  2  J 

Po  =  -T.  If) 


-sinP 


r2 


(13) 

(14) 

(15) 


where  K  =  K-k;  Ca  =  Ca/K^ . 

Similarly,  expressions  «j ,  Vj  and  pj  may  be  derived.  Since  these  expressions  are  very 
long,  they  are  not  reproduced  here. 

Using  equations  (13)  -  (15)  we  may  show  that  the  leading  order  evolution  equation  for  the 
film  rupture  is  given  by: 


^  ^  R± 

dx  2  ^  Fr^ 


sinP  +  /?e  • 


[-1 

L1-. 

d^h  - 

,"1 

iRe  ■ 

^  Ca 

3X3  ■"* 

\e/  J 

kAAT  ^*1  ^ 

’axj  ^BX 


Re^h^ 

3 


Ca ' 


[(£p2  +  e^(£p2) 


KAA^r 

B^h 

(A^+l)l 

^BhY\ 

+  /,iV+l[ 

9X2"^ 

h  ' 

.axJ  J 

(16) 


2  ax’ 


^smp  +  R. .  ^  -  1 )  A  ( (£V)  a  ^  )] 


kAN  dh 

1  ■  ax  ^ 
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subject  to  initial  conditions: 


h(X,0)  =  F(X) 


(17) 


Equations  (16)  and  (17)  may  be  solved  numerically  in  order  to  predict  the  stability  characteristics. 
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3.0  EFFECT  OF  ELECTROSTATIC  FIELD  ON  EVAPORATING  FILMS 

Here  we  consider  a  thin  liquid  film  joining  an  absorbed  film  and  the  meniscus  of  an  evano- 
ratmg  interface  under  the  influence  of  an  electrostatic  field.  From  the  previous  analysis  to  a  lead¬ 
ing  order  of  magnitude,  we  have  ^ 

A  J  +  (18) 

In  the  above  equation,  is  the  liquid  pressure;  the  vapor  phase  pressure;  A  the  Hamaker  con¬ 
stant  (negative  for  a  spreading  Uquid);  5  the  film  thickness;  the  fluid  dielectric  constant  and  E 
the  electric  field. 

Using  the  lubrication  approximation  for  the  liquid  flow  in  the  thin  film  we  have 


d^u  jPt 
~  dx 


with  boundary  conditions  given  by 


y  =  0:  M  =  0  and  y  =  5:  =  0: 


From  equations  (19)  and  (20)  we  may  write 

1  ^Px  ry2 

“  ~  a,  '  dx  ' 12 

The  mass  flow  rate  per  unit  width  is  given  by 


[f-a.] 


r  „  f  j  53  ‘'P. 

r  =  =  (Z 

0 

Following  Wayner  [47],  the  evaporative  flux  by  means  of  theKelvin-Claperon  equation  may  be 
written  3s 

”  =  <•  (r.v -T,)+b (p, -p„)  =  ■  la  (T,-T,)+b (p, -p,) ]  ■  (Z 

il+— 8j 
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ICO 


where 


=  evaporative  flux 


jr( 


\2nRT,X  RTJ,, 


\2kRT^XRT,, 


=  temperature  of  liquid-vapor  interface 
=  temperature  of  the  vapor 
A/i^  =  enthalpy  of  vaporization  per  unit  mass 

k  =  thermal  conductivity  of  the  liquid 

is  related  to  the  surface  temperature  by  the  one-dimensional  heat  conduction  equation 

(7^  _  )  =  mA/i^ .  Evaporative  flux  is  related  to  the  flow  rate  in  the  film  through  mass  bal¬ 

ance  in  the  following  form: 


=  -m 


From  equations  (23)  and  (24),  we  may  write 


1  d{  ,^Px\  1 


I  +  — 8 


where  AT  =  T^-T^ 


we  define  the  following  dimensionless  variables: 


53  =  -^ 
®0  aAT 


=  reference  thickness 


=  reference  pressure 


dimensionless  pressure  difference 
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=  scale  factor  for  X 


dimensionless  position 


dimensionless  thickness 


Dividing  equation  (18)  by  tCq  we  get  the  dimensionless  form: 


1  dr\^ 

=  (26) 


where  3  (^)  =  [(^) 

_1_  , 
■’'o'' 

(i-O 

where  eI  =  E^JF,  eJ  = 

EJ/F 

Equation  (25)  reduces  to 

1  flf. 

1 

"  (1+KTl)  ■ 

(27) 

Boundary  conditions  are: 

oo:  q  1; 

^  >  — oo 

(28) 

Equations  (26)  -  (28) 
rating  film  heat  transfer. 

may  be  solved  to  predict  the  effect  of  electrostatic 

field  on  the  evapo- 

I  = 


-A 
Iv^aAT 


s 

1  =  5: 


aAhmdfs 


K  = 


e  = 


(A) 
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4.0  EFFECT  OF  ELECTRICAL  FIELD  ON  EVAPORATING  FILMS 

In  the  previo^  two  cases,  it  has  been  assumed  that  the  coupling  of  the  electrostatic  problem 
to  the  fluid  dynamics  of  the  film  is  only  through  the  boundary  condition  imposed  at  the  liquid- 
yapor  interface.  Although  this  may  be  valid  in  the  case  of  weak  electrical  fields,  the  case  involv¬ 
ing  strong  electrical  fields  (which  is  of  practical  interest)  should  consider  the  proper  body  forces 
in  the  momentum  equation  and  the  work  terms  in  the  energy  equation. 

References  [11-20]  consider  the  effect  of  electrical  fields  on  single  phase  flows.  A  thorough 
literature  search  has  indicated  that  the  effect  of  electrical  fields  on  evaporating  films  or  two  phase 
flows  has  not  been  investigated  so  far.  ^ 

At  the  conclusion  of  the  present  task,  there  was  not  enough  time  to  formulate  this  problem 
in  detail.  Therefore,  only  the  governing  equations  for  the  problem  are  given  here. 

We  assume  that  the  fluid  is  incompressible.  The  current  in  the  fluid  is  assumed  to  obey 
Ohin  s  law  with  the  conductivity  a  function  of  temperature.  The  governing  equations  for  the  con¬ 
vective  heat  transfer  in  an  electric  field  may  be  written  as: 

Mass: 


V-  =  0 


(29) 


Momentum: 


.dr  ■ 


V- 


(30) 


where  Fvdw  is  the  volumetric  body  force  due  to  long  range  intermolecular  (Van  der  Wall)  forces. 
Energy: 


where 


dr  (t£2 

^  +  V.Vr=  aV27+^ 
or  pCp 


(31) 


p  =  density  of  the  fluid 


p.  =  viscosity 


q  =  charge  density 
E  =  electric  field 
V  =  velocity  field 
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E  “'’.-TPla^Jr 


=  liquid  pressure 


In  a  poorly  conducting  liquid,  the  currents  are  small  and  therefore  magnetic  fields  are  negligible. 
The  equations  governing  the  electric  field  are  given  by: 

q  =  V  -  (e£)  =  free  charge  density 
Vx£  =  0 


J= 


0 


J=<sE  +  qV 

where 

e  =  dielectric  constant 
♦ 

J  =  current  density 
<y  =  electrical  conductivity 

The  boundary  conditions  are  no  slip  at  the  solid  surface  and  the  same  as  previously  described  at 
the  liquid  vapor  interface. 
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